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Abstract 

Non-perturbative lattice studies of QCD in the chiral thermodynamic regime, where 
chiral symmetry is spontaneously broken, require to deal with almost quark zero modes 
in a theoretically clean and computationally efficient way. We discuss the basic features 
and some realistic tests of a formulation, known as lattice tmQCD, that fulfills these 
requirements. 
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1 INTRODUCTION 



Lattice field theory is known to provide a rigorous and systematically improvable way of 
studying QCD in the non-perturbative regime by means of Monte Carlo simulations. As 
in any numerical method, the coexistence of widely separated energy scales poses severe 
practical problems: for instance, the energy-momentum resolution a -1 , which is induced 
by the lattice regularization, should be kept much larger than both the typical hadron 
mass scale, say ~ 500 MeV, and the external momenta of the correlation functions that 
are computed. Moreover, the approximate and spontaneously broken flavour chiral 
symmetry that is exhibited by the strong interactions entails the need of performing at 
least a few simulations of lattice QCD in large volume and with small quark masses^]: 
a fully realistic setup would require a lattice of spatial volume L 3 , with M^L > 5 and 
M n ~ 140 MeV. 

The relevance of the approximate chiral symmetry for the low energy QCD am- 
plitudes can hardly be overestimated: it determines many aspects of the dynamics, as 
shown by the chiral effective Lagrangian approach, and strongly constrains the oper- 
ator mixings, which on the lattice may be particularly severe due to the dimensionful 
ultraviolet cutoff a -1 . The realization of the chiral symmetry on the lattice is known 
to be delicate since the pioneering work by Wilson [jij and in most cases the full flavour 
chiral symmetry is recovered only in the continuum limit. The highly remarkable ex- 
ception is represented by those lattice formulations where the Dirac operator satisfies 
the Ginsparg- Wilson relation which in turn entails an exact flavour chiral invari- 
ance at finite lattice spacing. The lattice regularizations with Ginsparg- Wilson quarks 
certainly simplify a lot the construction of renormalized operators, especially those rel- 
evant for the weak effective Hamiltonian, at the price however of a big overhead in the 
computational effort. In many cases one can avoid such an overhead by working within 
the framework of Wilson quarks, which seems to be more flexible and powerful than 
believed till few years ago, as we try to argue in the following. In general, the lattice 
regularization of the quark sector should be chosen with care depending on the physical 
applications and the related renormalization problems. 

In this contribution we focus on a lattice formulation of QCD, based on Wilson 
quarks and known as twisted mass QCD (tmQCD), that is particularly suited for dealing 
with the u and d quarks. After discussing in Section 2 how the simulations of lattice 
QCD account for the contribution of quark zero modes, we illustrate in Section 3 the 
basic features of lattice tmQCD. In Section 4 we briefly report on an exploratory non- 
perturbative study of 0(a) improved lattice tmQCD, which reaches a pseudoscalar to 
vector meson mass ratio of Mps/My — 0.47(1) and represents a successful test of several 
aspects of our approach. 



3 This task is significantly alleviated by studying the renormalization problems in finite volume by 
means of renormalization group and finite-size scaling techniques, e.g. in the Schrodinger functional 
scheme. 
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2 Lattice QCD and physical quark zero modes 

We assume that the reader is familiar with the lattice regularization of QCD introduced 
by Wilson Q and for simplicity we consider the theory with two mass-degenerate quark 
flavours. In this case the lattice action with Wilson quarks reads 

S = S g [U; gl] + a^xix) [(±7 • (V + V*) - a±V* • V + m' ) X ] (x) , (2.1) 

X 

where S g is the pure gauge action, x denotes the doublet of quark fields and v-w = v^w^, 
while = VjuffJ] and V* = V*[f7] stand for the forward and backward covariant 
derivatives on the lattice. The hard breaking of the flavour non-singlet axial generators 
induced by the Wilson term can be compensated up to 0(a) effects by tuning the 
parameter m' and the coefficients that parameterize any chirality-violating operator 
mixings ||. The leading cutoff effects can be removed via the on-shell 0(a) improvement 
Q a la Symanzik. 

2.1 Valence and sea quarks 

In any correlation function to be evaluated via Monte Carlo simulations, due to the 
huge dimensionality of the vector space spanned by the Grassmann variables in the 
Euclidean path integral, it is customary to integrate out analytically the fermionic 
degrees of freedom. In the case of the correlator of two local pseudoscalar densities with 
isospin index a = 1, 

C&(x -y) = -(P l {x)P\y)) = -Z- 1 j dt/d X dxexp(-S) [xi^x]{x)[xi^x](y) , 

(2.2) 

this procedure yields 

CH(x-y) = J dU P[U;glm' q ] ■ 

■ tv{[D W)C + m' q ]~ 1 (x,y)'y 5 ^-[D WjC + m q ] -1 (y, xhs^} , 
P[U;glm' q } = Z' 1 exp(-S g [U, g 2 }) det({Ay,c + m' q }[U}) > , (2.3) 

wheref] Z denotes the Euclidean partition function and -Dw,c = %, c [f/] is the critical 
(two-flavour) Wilson-Dirac operator: 

Av,c = h-( V + V *)-f V *- V + m c (2.4) 

The parameter m' q = m' — m c is hence proportional to the renormalized quark mass 
that appears in the PCAC Ward identity: the massless theory is obtained for m' = 
m c (gQ). Gauge configurations U are generated via suitable algorithms with probability 
P[U; gQ,m' q ], and on each configuration U the quark propagator from the lattice site x 
to the site y, [D\v, c + mq] _1 (y,x), can be computed -for fixed x- by solving a linear 
system. By expanding the correlator (|2.2| ) around the trivial perturbative vacuum, one 

4 The symbol tr{. . .} denotes the trace over flavour, colour and spin indices, whereas the symbol 
det{. . .} stands for the determinant with respect to all indices, including the space-time ones. 
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can easily check that in eq. (|2,3| ) the trace term involving quark propagators corresponds 
to valence quark diagrams dressed with any kind of purely gluonic corrections, whereas 
the term det(Z?w,c + m 'q) accounts for the sea quark corrections to the aforementioned 
dressed valence quark diagrams. 

2.2 Quenched and partially quenched lattice QCD 

It is well-known that with the established simulation techniques the simulation of the full 
theory, including the sea quark effects, has a very high computational cost, which quickly 
increases as the pion mass is decreased towards realistic values. On the other hand, it 
is technically trivial to choose different values for the parameter m' q that appears in the 
fermionic determinant and its counterpart in the inverse Dirac operator: 7^ m' al . 

A moment of thought reveals that such a modification of lattice QCD corresponds to 
a statistical model with extra spin-1/2 ghost fields, which is local and renormalizable 
by power counting but violates reflection positivity (see e.g. Ref. ||). The general case 
m qvai m q,sea * s referred to as partially quenched QCD, whereas the particular case 
l m qvail < l m q,scal = 00 corresponds to the well-known quenched approximation. 

It should be noted that these approximations in general break down as the valence 
quark mass Tnf va [ — ► in the thermodynamic limit. In this regime the flavour chiral 
symmetry is spontaneously broken: if the chiral condensate is to be non- vanishing, 
the gauge configurations carrying zero modes of the Dirac operator must receive a 
finite weight in the Euclidean path integral, even in the limit (taken at infinite spatial 
volume) mq Sea = va j ► Q . It is clear that in the full theory the integration over 
the fermionic variables on any gauge background can not yield divergences: in presence 
of quark zero modes the infinities in the quark propagators must be compensated by 
the zeros in the fermion determinant, see eq. Q2.3| ). Such a delicate compensation is 
no longer guaranteed if the condition m' q sea = m' q val is violated: fermionic observables 
may hence diverge on gauge configurations carrying quark zero modes, which causes the 
breakdown of the quenched and partially quenched approximations. 

However, when working sufficiently away from the thermodynamic chiral limit, the 
quenched -or partially quenched- approximation is expected to be reasonably accurate, 
at least for those quantities that are not very sensitive to sea quark effects. An example is 
given by the ratios of hadron masses J7[|, with the ?/-meson mass being the most striking 
exception. Indeed, given the high computational cost of simulating the full theory, the 
quenched approximation has been widely used as a testbed for lattice techniques and 
for first non-perturbative estimates of quantities such as renormalized couplings, hadron 
masses and matrix elements, order parameters of phase transitions. The chiral effective 
Lagrangian for quenched QCD has also been worked out M with the aim of identifying 
and parameterising the deviations from the full theory close to the chiral limit. 

Analogous remarks apply for the partially quenched approximation, the quality 
of which depends on the ratio m' qfiea / m' aV In particular, it has been remarked that 
the low-energy (Gasser-Leutwyler) constants of partially quenched QCD with Af quark 
flavours coincide with those of the fully unquenched theory, provided that all the quark 
flavours are light enough for the chiral perturbation theory to be applicable ||. In 
the physically relevant cases, which are Af = 3 and -to some extent- Af = 2, varying 
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m qvai w hile keeping fixed m qsea can then be very convenient, since it allows to investi- 
gate the dependence of the observables on m! val without performing many unquenched 
simulations. 



3 Lattice tmQCD 

Lattice tmQCD is an extension of the widely used formulation with Wilson quarks, 
from which it differs in that the physically non-vanishing quark mass term is in general 
not aligned with the Wilson term in the flavour chiral space. This simple modification 
brings definite advantages concerning the simulations with light quarks -especially in 
the quenched or partially quenched case- and the renormalization properties of some 
phenomenologically important quantities, such as the leptonic decay amplitude of pseu- 
doscalar mesons or the mixing amplitude in the K°-K system PJlO|,ll|. 



3.1 Lattice Wilson quarks and spurious zero modes 

While the considerations of Subsections 2.1-2.2, apply to any sensible lattice regular- 
ization of QCD, working with Wilson fermions entails a further technical problem that 
renders particularly difficult -or even impossible- the simulations with light quarks. The 
problem arises whenever the quenched (or partially quenched) sample of configurations, 
as determined by a given choice of the bare parameters, includes gauge backgrounds on 
which the critical Wilson-Dirac operator -D\y,o eq. (2.4), has one or more eigenvalues 



with negative real part Re(A) < and (almost) zero imaginary part. Under these con- 
ditions, the massive Dirac operator -Dw.c + m q va i i s singular for values of m' qvaJ > as 
soon as "i qva i + Re(A) = for some of the real negative eigenvalues A. 

The fermionic observables, which involve [-£>w,c + m qv ai] _1 ' ma y receive (almost) 
divergent contributions on the aforementioned gauge backgrounds, spoiling the expected 



decrease of the statistical errors with the number of independent measurements [12]. 
As an example of this phenomenon, we show in Fig. |l| the Monte Carlo history of the 
normalized relative standard deviation for the pion channel correlator / = /p (xo) at 
xq = 24a = T/2, see Section 4.1. The normalization of the standard deviation is such 
that it should approach a constant in the limit of infinite statistics: in the case of 
the standard Wilson regularization, which corresponds to the open symbols in Fig. |], 
the problem is apparent. Via the combined effect of lattice artifacts and statistical 
fluctuations, "spurious" quark zero modes^] anticipate at non- vanishing values of m' q 
the breakdown of the quenched or partially quenched approximation. 

In practice, in the quenched case the rate of occurrence of gauge backgrounds with 
spurious quark zero modes -also called "exceptional configurations"- depends on the 
values of m q val and g$ , on the physical linear size L of the lattice and on various details 
of the lattice regularization. The rate tends to increase when decreasing m qval and 
increasing g$ and L/a, as well as when switching on the coefficient, c sw (5o)> of the coun- 
terterm that is needed for the on-shell 0(a) improvement of the fermionic action. For 
instance, in the regularization with plaquette gauge action and non-perturbatively 0(a) 

5 These spurious quark zero modes should not be confused with the quark zero modes that play an 
important physical role in the thermodynamic chiral limit of renormalized QCD. 
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improved Wilson quarks, the problem is felt for L > 1.5 fm at values of corresponding 
to a ~ 0.1 fm and at valence quark masses that are about half the strange quark mass. 
A similar problem is also expected in partially quenched simulations, although with a 
lower rate depending on the ratio m q va i/ m' q ^ sea , an d has indeed been observed, see e.g. 



Ref. |J. 



In the fully unquenched case, one can expect troubles at the algorithmic level 
with light Wilson quarks on rather coarse lattices. This is because the state-of-the- 
art algorithms implement stochastically the fermion determinant that appears in the 
probability measure for the gauge configurations: almost exceptional configurations may 
hence be proposed, but are then almost certainly rejected. In this process however the 
simulation algorithm, e.g. the standard HMC one^|, undergoes a severe slowing-down, 
due to a decrease in the acceptance rate and an increment of the condition number of 
the Dirac matrix before the accept/reject test. 

3.2 Action and symmetries 



As already known since 1989 [15], the problem with spurious quark zero modes is absent 



in the two-flavour theory if one considers the action 

S w [U,ip,ip] = 5 g + o 4 ^] ^(x)[(I>w,c + m q + j/J q 75i" 3 )^](s) , (3.1) 

X 

where tp is a flavour quark doublet, the matrix r 3 acts in the flavour space and the 
boundary conditions for finite- volume systems may remain unspecified for a while. Since 
m c is known in practice with finite precision, the exactly known bare parameters are g$ , 
/i q and mo = m c +m q . It is easy to see that, as long as /x q ^ 0, on any gauge background 
the lowest eigenvalue of the Hermitean square of the matrix (-D\y,c + m q + i/i q 75T 3 ) is 
bounded from below by (a/x q ) 2 . 

The action ( |3.1| ) represents a sensible regularization of QCD with Nf = 2 mass- 
degenerate quark flavours, but in a quark field basis that is chirally twisted with respect 
to the standard one |J. To illustrate this point, let us focus on the simple case m q = 0, 
Hq ^ 0, and consider the very same lattice theory in the standard quark field basis, 
which is obtained by a suitable axial rotation with generator r 3 : 



<t-3 \ / ~-3 \ // 7T 

X = exp ( iuyys-] ip , x = 4> exp \ , w = arctan— 2- = - 



(3.2) 

The action ([O]) then reads 

Sw[U,x,x] = 5 g + fl 4 ^xW[(A m w,c +/j,q)x](x) , (3.3) 

x 

where 

Amw.c = h ■ ( V + V *) + ^t75T 3 V* • V - i 75 r 3 m c . (3.4) 

The connection, eq. (|3.2| ), between the two lattice quark bases makes obvious that 
the chiral limit is obtained for m q = /i q = with m c being the usual function of 

6 Evidence for a large increase of the fermionic force at small quark mass values is reported in Ref. |F3 . 
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og. Inspection of the symmetries of the action (^3) shows that in the chiral limit one 
vector and two axial generators -out of the six generators of the flavour chiral group- 
are preserved by the lattice regularization, while parity is preserved only up to a flavour 
exchange (Pp symmetry) and all other symmetries are as usual with Wilson fermions. It 
follows that /x q is multiplicatively renormalized, /iR = Z^g^) /i q , while = Z^gfy m q . 
Power counting renormalizability and the recovery of the full flavour chiral symmetry in 
the continuum limit Q, together with the exact Pf invariance, imply that parity must 
also be recovered in the continuum limit. In contrast with other approaches based on 



the action ( |3.1[) [12], the correlation functions of the massive QCD are obtained at finite 
/i q , so that the problem with the spurious quark zero modes is certainly solved. 

In the general case, where both m q and /z q are non-vanishing, the situation is 
fully analogous, but the relation between the two lattice quark bases involves an angle 
u 7^ it/2. However, owing to the complications arising from the Wilson term, it is 
convenient to perform first the renormalization (and possibly the 0(a) improvement) 
of the correlation functions in a given quark basis and then transform to the "physical" 
basis, i.e. the one where the quark mass is coupled to the singlet scalar density. The 
first step implies that the continuum limit is approached at fixed values of <? R , thr, /xr 
and the normalization conditions for the composite fields. Concerning the second step, 
for a wide class of renormalization schemes, the relation between the renormalized fields 
in two different bases that are related via a non-singlet axial rotation takes the same 
form as at the classical level ||. In terms of polar quark mass coordinates, 



tan a = — ^ , Mr = J ok + m£ , (3.5) 
tor V n 

the angle a identifies the quark basis and hence specifies the axial rotation with generator 
r 3 /2 that allows to combine the renormalized correlation functions so to obtain finite 
correlators with given continuum quantum numbers]]. On the latter correlators the 
partial breaking of isospin and parity is an effect of order a/iR and represents no serious 
drawback. 

If one identifies the two quark flavours with the u and d quark, the tiny mass differ- 
ence between them can safely be taken into account by means of chiral perturbation the- 
ory. Alternatively, lattice tmQCD can also be formulated for a doublet of non-degenerate 



quarks, whilst retaining the protection against spurious quark zero modes [17]. Heavier 
flavours of Wilson quarks can be added e.g. in the usual way to the lattice tmQCD 
action for Nf = 2. 

4 A test-study of lattice tmQCD with light quarks and in large volume 

We present here some preliminary results of a high-statistics exploratory study of lattice 
tmQCD in the thermodynamic chiral regime. The study aimed at testing the absence of 
exceptional configurations, the computational cost and the magnitude of cutoff effects 
for a few typical observables. We adopt in the following the notation of Ref. [19|, to the 



equations of which we refer with the prefix " I" , and postpone many technical details to 
a forthcoming publication [2fl]. 



7 This point of view might lead to a reinterpretation of the "spontaneous breakdown of parit y an d 
isospin" in lattice QCD with Wilson fermions 



16j, which was observed by employing the action (3.1) 
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4.1 Observables and simulations 



In this test-study we choose to work in the quark basis of action ( |3.1[ ) and implement the 
non-perturbative 0(a) improvement of the action and the relevant operators along the 
lines of Ref. [|l8|. Attention is restricted to the pseudoscalar and vector meson masses, 
Mps and My, the pseudoscalar (leptonic) decay constant, Fps, and the polar quark 
masses, Mp and a, which are defined according to eqs. (I.3.3)-(I.3.20) for /ur, and mR 
and eq. (|3~5|). 



Set, (i 


L/a, T/L 


Apneas 


m R /M R 


M R r 


M PS r 


F PS r 


M v r 


Al, 6.0 


16, 2 


650 


-0.016(3) 


0.2729(15) 


1.711(7) 


0.455(5) 


2.662(40) 


Al',6.0 


16, 3 


650 


-0.016(3) 


0.2729(15) 


1.714(6) 


0.455(6) 


2.656(42) 


A2, 6.2 


24, 2 


535 


-0.014(2) 


0.2558(16) 


1.623(8) 


0.456(5) 


2.557(32) 


B1L, 6.0 


24, 2 


260 


0.017(3) 


0.1949(11) 


1.452(6) 


0.432(6) 


2.517(35) 


Bl, 6.0 


16, 2 


535 


0.001(3) 


0.1949(11) 


1.455(8) 


0.428(5) 


2.513(47) 


B2, 6.2 


24, 2 


300 


-0.004(4) 


0.1962(12) 


1.420(9) 


0.436(7) 


2.462(41) 


C, 6.0 


24, 2 


260 


0.083(5) 


0.1205(7) 


1.160(6) 


0.401(6) 


2.485(59) 



Table 1: Statistics and renormalized quantities obtained in our simulations, which are 
identified by a label and the value of (3 = 6/oq. The statistics is specified by the number 
of measurements, A meas , on almost independent gauge configurations, while the values 
of mo and fi q can be found in Ref. [20]. 



We work in the quenched approximation considering systems of physical size L 3 T 
with L such that MpgL > 4.5 to suppress finite volume effects. In practice we take 



L = 1.5 fm or L = 2.2 fm and T/L = 2 -r 3. Following Ref. [19], we impose boundary 



conditions of Schrodinger functional (SF) type and compute the SF correlators 

/AVo), ti\x ), fy 2 (x ), *#(*«,), kf{x Q ), ft 1 . (4.1) 

We hence construct the linear combinations of renormalized and O(a) improved SF 
correlators that correspond to operator insertions (at time xq) with well defined con- 
tinuum quantum numbers, see eqs. (1.3.10)— (1.3.11) with a given by eq. ([3.5|). Namely, 
the correlators [/A/]r(^o) an d [fp}]n(xo) correspond to the insertion of the isotriplet 
pseudoscalar operators (^4p)o and (-PpJ 1 , while the correlators [/cy 1 / ] R (xo) and [fc^] R (xo) 
correspond to the insertion of the isotriplet vector operators (V^)^ and (IpJjL. In the 
limit of large xo and large (T — xq) and up to cutoff effects, these non-vanishing cor- 
relators are expected to be dominated by the lowest isotriplet pseudoscalar and vector 
meson states. 

An overview of our simulation parameters, statistics and preliminary results for 
renormalized quantities is given in Table 1. The renormalized gauge coupling g 2 ^ is 



eliminated in favour of the length scale ro [21], which is known to be about 0.5 fm, 
while the lattice spacing value corresponding to = 6 (6.2) is a ~ 0.093 (0.068) fm. 
Our most critical simulation (set C, (3 = 6), where we employed a CGNE solver for 
the SSOR-preconditioned version of the Dirac matrix (Z?w,c + "iq + ifJ-ql5i~ 3 ), required 
~ 230 GFlops x day. 
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50 100 
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Figure 1: The square root of the relative a priori variance of / = fp (xo = 24a) 
versus the number of measurements N: data from the simulation C (filled symbols) and 
another simulation with the same values of (3 and Mr, but a = (open symbols). 



4.2 Results 

We find that lattice tmQCD allows, as expected, to safely work in a region of parameters 
which would be inaccessible with ordinary Wilson quarks: see e.g. Fig. |l| as well as the 
findings of Ref. [24]. For a given number of independent measurements, the statistical 



errors on Mps and My are comparable, up to a factor of one to three, to those found 



e.g. with domain wall quarks [22]. The CPU time effort, e.g. for the data sets Al, Al' 
and A2, is in line with the computational cost for ordinary Wilson quarks. 

The partial breaking of parity and isospin that is peculiar of lattice tmQCD is found 
to be a minor problem within our small statistical errors. In this respect it should be 
noted that we work at small values of a/i q , namely 0.0266 > afx q > 0.0117, and consider 
observables in physical channels where the lowest state is lighter than the lowest state of 
the corresponding channels with flipped parity and isospin. While deferring the details 



of our analysis to Ref. [20|, we show in Fig. ^ an example of effective masses extracted 
from SF correlators, where the correlator [/^]r(xo) receives contributions of order a/i q 
that are visible at large xq. Analogous effects are expected and found to be negligible 
within statistical errors for both [/pf]R(xo) and the vector channel correlators. 

As detailed in Refs. |18,10|, we expect the relations among our observables and the 
renormalized parameters ro and Mr to be O(a) improved. In particular, when working 
at a = 7r/2 + O(a), which is the case of our study, the knowledge of a few counterterms 
(those with coefficients Z g , m c and c sw ) suffices to obtain an O(a) improved estimate of 
-Fps- In order to check for the residual scaling violations, we produced data at (3 = 6.2 
(sets A2 and B2), while keeping a, Mr and ro fixed. The small mismatch in Mrto for 
the set A2 was corrected by employing estimates of the dependence of our observables 
on M R r . 



We also reanalysed the data of Ref. |23|], which were produced at (3 = 6, 6.1, 6.2, 6.45, 
by imposing precisely the same renormalization conditions as in this study of tmQCD. 
We then performed a continuum extrapolation of these data, assuming a purely quadratic 



S 





0.29 




0.27 




0.25 






o 


0.23 


5 




i 


0.21 




0.19 




0.17 




0.15 







plateau window 
plateau window 



1 1 



10 15 



in; 



or [f A ,] R 
or [f p ] R 



115 



20 25 

xja 



30 35 40 45 



Figure 2: Pseudoscalar effective masses extracted from the correlators [/p]r = [/p/]R 
and [/a']r = [/a']r f° r the simulation C. The tiny circles denote our (good) fit to the 
effective masses from [/a']r: at large xq a peculiar 0(a/x q ) contribution is present. 



dependence on (a/ro) 2 and discarding the data at /3 = 6. However, the resulting es- 
timate of Fps is not fully 0(a) improved, as for one of the necessary improvement 
coefficients, 6a(<7o)> orn y the one-loop estimate could be used. The outcome of this ex- 
ercise is compared with the results from tmQCD in Figs. 3-4, omitting the case of Myro 
where cutoff effects are hardly visible within statistical errors. The estimators of Mps 
and Fpg that are obtained from lattice tmQCD show rather small cutoff effects, which 
agrees with the findings of a scaling test [19] in intermediate volume (L = 0.75 fm). 



5 Conclusions 



Lattice tmQCD is well suited to perform non-perturbative studies of QCD in the chiral 
(thermodynamic or finite- volume) regime for all cases where it is technically sufficient to 
recover chiral symmetry in the continuum limit. The framework has been successfully 
tested in the quenched approximation and can straightforwardly be extended beyond 
it. 

The tmQCD project is part of the ALPHA Collaboration research programme. We 
acknowledge the very pleasant collaboration with P.A. Grassi, S. Sint and P. Weisz, as 
well as fruitful discussions with M. Liischer, G.C. Rossi and R. Sommer. 
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